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We present a dynamic model based on the linear sigma model with constituent quarks 
that allows for studying the explicit propagation of fluctuations at a chiral critical point and 
a first-order phase transition. The coupled dynamics of the sigma field and the quarks is 
derived selfconsistently. Hereby, the sigma field evolves according to a semiclassical Langevin 
equation of motion, while the evolution of the quarks is described fluid dynamically in order 
to model the expansion of a hot fireball in heavy-ion collisions. Dissipative processes and 
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fluctuations in the Langevin equation are allowed under the assumption of energy-momentum 

conservation of the coupled system. We study the evolution of the constituent quark masses 

C*l i in this nonequilibrium set up. 
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I. INTRODUCTION 

The QCD phase diagram is expected to exhibit a rich phase structure. At larger baryochemical 
potential, as achieved at the upcoming ¥K. 



transition is expected from model studies 
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R project [l| at GSI Darmstadt, a first order phase 
-J]. Interesting observables could here be based on the 
growth of fluctuations due to the nonequilibrium effect of supercooling leading to nucleation and 
spinodal decomposition [5-8]. 

At zero baryochemical potential the nature of the phase transition of QCD is well understood 
from lattice QCD calculations, which show that it is an analytic crossover [91]. As a consequence 
there must be a critical point, which terminates the line of first order phase transitions. Unfortu- 



nately, there are no analytic techniques to solve the QCD partition function in this regime. Thus, 
one has to rely on model calculations and explore the phase diagram in experiment. In equilibrium 
systems fluctuations and correlations of the order parameter diverge at the critical point. Cou- 
pling particles to the sigma field, the order parameter of chiral symmetry, leads to a nonmonotonic 
behaviour in fluctuations of net-charge or net-baryon number multiplicities 
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111 ] . The key ingre- 



dient is the correlation length which becomes infinite in a system at a critical point. In a realistic 
evolution of a heavy-ion collision, however, the growth of the correlation length is limited by the 
size of the system and by the finite time, which the dynamic systems spends at a critical point. 
Relaxation times also become infinite at the critical point, a phenomenon called critical slowing 
down. Even if the system is in equilibrium above the critical point it is necessarily driven out of 
equilibrium by passing trough the critical point. Assuming a phenomenological time evolution of 
the correlation length with parameters from the 3d Ising universality class it was found that the 
correlation length does not grow beyond 2 — 3 f m [12]. 

The explicit propagation of fluctuations coupled to a dynamic model is a necessary step to- 
wards understanding the QCD phase diagram from heavy-ion collision experiments. In chiral fluid 
dynamic models [131 . Il4j the propagation of the fields in the chiral sector is coupled to a fluid 
dynamic propagation of the constituent quarks. The expansion and cooling of the fluid does thus 
drive the underlying model through the phase transition. In the following we present the model 



of nonequilibrium chiral fluid dynamics 15j] with a focus o the evolution of the constituent quark 
masses. 

II. NONEQUILIBRIUM CHIRAL FLUID DYNAMICS 

Starting from the quark meson model [16] the coupled dynamics of the order parameter of chiral 
symmetry, the sigma field, and the fluid dynamic expansion of the quarks have been derived [la ]. 
The Langevin equation for the sigma mean-field reads 

d^a+ 8 -^+r,d t a = i. (1) 

da 

The effective potential to one-loop level is given by 



V eS (a,T) = U(a)-2d q Tj^]n(l + e^(-^ , 



(2) 



The classical potential U in ([2]) is of Mexian hat shape 

A 2 



Via) =^{a 2 -u 2 ) 2 -h q a-U . (3) 
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FIG. 1: The equilibrium quark masses as a function of temperature for a critical point and a first order 
phase transition 



The parameters are chosen such that chiral symmetry is spontaneously broken in the vacuum, 
where (a) = f n = 93 MeV. Chiral symmetry is broken explicitly by the term h q a. In the effective 
potential V e s the quark energy is given by 



E 2 = p 2 + ml 



(4) 



and the quark mass is generated dynamically at the phase transition, where the sigma develops a 
finite expectation value and thus m q = g(cr). As due to the explicit symmetry breaking term (a) 
does not vanish exactly in the chirally restored phase, but has a small finite value, m q does not 
exactly vanish either. The effective potential is given for fis = 0, where the correct vacuum value 
for the constituent quark mass is obtained by a coupling of g = 3.3 between the sigma field and the 
quark fields. For a qualitative study we fix the baryochemical potential at /i# = and tune the 
phase transition by changing the coupling constant. For g = 3.63 one finds a corresponding critical 
point at T c = 139.88 MeV by the vanishing curvature of the effective potential at the minimum and 
a first order phase transition for g = 5.5 and T c = 123.27 MeV by the appearance of two degenerate 
minima. For these couplings the phenomenologically known value of the constituent quark mass 
comes out wrong. We, however, focus on the qualitative behaviour and leave the extension of the 
model to finite \jlb to future work. In figure HT1 we. therefore, show the equilibrium quark masses 
devided by the vacuum quark masses. For a first order phase transition we see a large discontinuity 
at the phase transition due to the two degenerate minima, while the transition is smooth at the 



critical point. The damping term is 
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FIG. 2: Evolution of the constituent quark masses for a scenario with a critical point and a first order phase 
transition (left). The temperature evolution is shown for comparison (right). 
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V = \ 2.2/fm for 2m q > m a > 2m. K 

for m a < 2m n ,2m q 

The stochastic field in the Langevin equation ([1]) has a vanishing expectation value 

<£W> ? = o, 

and the noise correlation is given by the dissipation-fluctuation theorem 

mW)) ^ = ^(t ~ t')m aV coth {^) . 
The local pressure and energy density of the quarks is given by 

p(a, T) = -V cS (a, T) + U(a) e(a, T) = T^^l - p(a, T) . 



(5) 



(6) 



(7) 



(8) 



In the relativistic fluid dynamic equations we find a source term S v allowing for the energy dissi- 
pation from the system to the heat bath 



0^ u = S u . 



(9) 



III. DYNAMIC QUARK MASSES 



We solve equations ([T]) and Q numerically for simple initial conditions: an almond-shape, 
boost-invariant distribution of the energy density in the transverse plane and longitudinal direction 
respectively. First numerical results have been presented in [20|. In figure HU (left) we present the 
time evolution of the constituent quarks mass during the expansion of the quark fluid and the 



relaxation of the sigma field for both scenarios, with a critical point and with a first order phase 
transition. For a better understanding we also show the average temperature of the quark fluid in 
the same rregion. The transition temperature is crossed around t ~ 5 fm for a first order phase 
transition and around t ~ 4 fm for a critical point. At the first order phase transition we observe 
the interesting effect of reheating between 5 and 10 fm. In [2l| we discussed how this effect is 
caused due to the explicit energy-momentum conservation of the coupled system. We see that 
during the evolution of the system the constituent quark masses increase to their vacuum value. 
For a critical point scenario this behaviour is rather smooth and begins even at temperatures 
above the phase transition as the effective potential changes its shape gradually. During a first 
order phase transition, however, the high-temperature expectation value is separated from the 
vacuum expectation value by a finite barrier. Parts of the system are trapped in the unstable 
minimum below the phase transition. The jump of the constituent quark masses at the first order 
phase transition are thus smirred out and deviations from the equilibrium form are larger than in 
a critical point scenario. 

IV. SUMMARY 

We have presented a model of chiral fluid dynamics including dissipation and noise based on 
the linear sigma model with constituent quarks. The latter receive their vacuum mass by a finite 
expectation value of the sigma field in the chirally broken phase. We and studied the time evolution 
of the constituent quark mass. Due to supercooling at the first order phase transition we observe 
a retarded approach to the vacuum value. There are substantial parts below the phase transition 
where the quarks in the fluid have small masses compared to mj =0 . It will be interesting to study 
the influence of this effect in coupling to the confinement-deconfinement phase transition. Work 
including the Polyakov-loop is in progress 22]. 
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